#ifndef _MCMC_LOW_LEVEL_STAT_4Affinity_H
#define _MCMC_LOW_LEVEL_STAT_4Affinity_H

#include "mcRaster.h"
#include "mcmcHistogram.h"

class MCMC_FilterResp;

class MCMC_LowLevelStat4Region
{
public:
	// construction & destruction
	MCMC_LowLevelStat4Region();
	~MCMC_LowLevelStat4Region();
	void Free(void);

	MCMC_LowLevelStat4Region & operator=(const MCMC_LowLevelStat4Region &a);

	// operations
	void	SetNumFilters(const int num);
	void	Clear(void);

	void	Update(MCMC_FilterResp &resp, MCMC_PixelSet &pixel_set);

public:

	Histogram		*p_histogram_obs;
	Histogram		his_intensity;

	double	mean;
	double	var;
	int	num_pixel;

	double	a,b,c;	// fit it with a plane

protected:
	int		num_his;
};





#endif
